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Abstract — The interference channel achievable rate region is 
presented when the interference is treated as noise. The formula- 
-tion starts with the 2 user channel, and then extends the results 
■to the n user case. The rate region is found to be the convex hull 
' of the union of n power control rate regions, where each power 
control rate region is upperbounded by a (n — 1) -dimensional 
hyper-surface characterized by having one of the transmitters 
transmitting at full power. The convex hull operation lends itself to 
a time-sharing operation depending on the convexity behavior of 
those hyper-surfaces. In order to know when to use time-sharing 
rather than power control, the paper studies the hyper-surfaces 
convexity behavior in details for the 2 user channel with specific 
results pertaining to the symmetric channel. It is observed that 
most of the achievable rate region can be covered by using simple 
On/Off binary power control in conjunction with time-sharing. 
The binary power control creates several corner points in the 
n dimensional space. The crystallized rate region, named after 
its resulting crystal shape, is hence presented as the time-sharing 
convex hull imposed onto those corner points; thereby offering a 
viable new perspective of looking at the achievable rate region of 
the interference channel. 



I. Introduction 

One important communication model in wireless communi- 
cation is the interference channel, which is subject to intensive 

■ research nowadays. For example, the model is relevant for 
'cellular networks in which multiple base stations transmit data 
• to their respective subscribers and thereby causing interference 
[at the unintended receivers, and ad-hoc networks in which 

■ nodes are active at the same moment in the same frequency 
band. For a better understanding of the interference channel, 
it is crucial to know its capacity region, i.e., the maximum 
I set of all achievable rate points. It serves also as a benchmark 
for the comparison of different schemes. Unfortunately, the 
capacity region of the 2— user interference channel has been 
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an open problem for about 30 years jS), ||4l. Information- 
theoretic bounds through achievable rate regions have been 
proposed, most famously with the Han-Kobayashi region Q. 
The capacity of the Gaussian interference channel under strong 
interference has been found in (|6l, l?). Recent results on the 
2— user interference channel to within one bit of capacity have 
been shown in fBl, where a simplified Han-Kobayashi scheme 
was used in which the message is split in two parts, the 
private part and the common part. The transmit signal is then 
a superposition of those two signals. By a smart allocation of 
power between those two parts it was shown that this scheme 
is asymptotically optimal using a new metric, which is referred 
to as the generalized degrees of freedom. The generalization 
of the obtained results to the n-user case is rather difficult. As 
such, for the n-user case mainly the capacity slope as a function 
of the Signal-to-Noise-Ratio (SNR) are known. It was shown 
in t9j that for very high SNR, the capacity can be approximated 
by C = I \og{SNR) + o {\og{SNR)), where the second term 
vanishes by definition for extremely high SNR. 

The aforementioned referenced literature focused mainly on 
the 2— user interference channel from an information-theoretic 
point of view with highly sophisticated and thus quite com- 
plicated transmitters and receivers. There are other works in 
literature that tackle the practical issues in order to improve 
the performance of the interference channels. Power control 
is one element of critical importance. In IfTOl . a framework 
for the uplink power control is constructed and iterative power 
control is proposed. Adaptive modulation and coding (AMC) 
can be combined with power control to enhance the network 
performance ifTTI . For multiple channel (such as OFDM) and 
multiple cell case, joint AMC and power control have been 
widely employed lfT2l . Beamforming and spatial diversity can 
also be utilized when communicating over the MIMO channel 
ifTSl . Interference avoidance lfT4l has also attracted many recent 
attentions. Finally, many distributed solutions are proposed 
ifTSl - lfTTl with the benefit of simple implementation or low 
data overhead. 

In this paper, the achievable rate region is discussed for the 
?7 — user interference channel when the interference is treated as 
additive Gaussian white noise and no multi-user detection is 
employed. Examples where we encounter the need to define 
such rate region are found in multicell communications, in 
addition to mesh and sensor networks where the preference 
is to use low-complexity transceivers. It is also interesting 
to note that using the strategy of treating interference as 
additional noise proves to be asymptotically optimal, i.e., it 
was shown in lH that treating interference as noise is optimal, 
as long as the interference power in dB is lower than half 
of the useful signal power given that the power to noise 
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Fig. 1: n— user interference channel 



ratio is asymptotically high. This result was extended to the 
non-asymptotic case independently by three research groups 
lfT8l - ll20l . The generalization to the asymmetric case and the 
n— user interference channel is given in lfT8l - ll20l as well, 
where it was shown that it is optimal to treat the interference 
as noise whenever a similar (sufficient) condition holds. The 
references |fT8l - ||20l considered only the Gaussian interference 
channel, while the general discrete many-to-one and one-to- 
many memoryless channels were investigated in 1211 . It was 
shown in |21 1 that treating interference as noise is also optimal 
in the discrete memoryless channel as long as the received 
signal at the interfered receiver is stochastically degraded 
compared to the received signals of the other receivers. The 
optimality of treating interference as noise for the multiple 
antenna case has been considered by Il22l - ll24l . 

This paper finds the achievable rate region for the 72— user 
interference channel as the convex hull of the union of n 
rate regions formed via power control, where each rate region 
is upperbounded by a hyper-surface of dimension n — 1 
characterized by having one of the transmitters operating at 
full power. Given that there is a convex hull operation imposed 
onto the hyper-surfaces, it is important to know their convexity 
behavior in order to determine when time-sharing should be 
applied. This is treated in details for the 2— user interference 
channel. As the convex hull operation lends itself naturally to a 
time-sharing operation, and based on the convexity conditions 
found, the paper discusses when a time-sharing strategy should 
be employed rather than pure power control, and then presents 
specific results pertaining to the 2— user symmetric channel. It 
is observed that the achievable rate region can be practically 
approximated by using simple On/Off binary power control 
in conjunction with time-sharing. The On/Off binary power 
control creates several corner points in the ri— dimensional rate 
region, and employing a convex hull time-sharing operation 
on those points achieves what is denoted as a crystallized rate 
region. 

The system setup is presented in section |ll] Section |lll] 
discusses the achievable rate region for the 2— user interference 
channel, and then generalizes the results to the user case. 
Section |IV] focuses on characterizing the 2— user rate region 
in terms of convexity or concavity and when time-sharing is 
optimal with specific results to the symmetric channel. Section 
rvl introduces the concept of the crystallized rate region where 
time-sharing and On/Off binary power control are used. Finally, 
the conclusion is drawn in section [Vl] 



II. System Model 

The n— user interference channel is presented in Fig[T]with 
n transmitters and n receivers. The i*'' transmitter transmits 
its signal Xi to the intended z*'* receiver with power Pj. The 
receivers have independent additive complex white Gaussian 
noise with zero mean and variance of cr^j. Each transmitter 
is assumed to have a peak power constraint of Pmax- Each 
transmitter has a single antenna and communicates over a 
frequency flat channel, gij denotes the channel power gain 
received at the i*'* receiver from the j*'' transmitter Therefore, 
(7i i is the channel gain of the i*'^ desired signal, whereas 
gij with j i represents the interfering channel gain. P is 
the transmit power vector of length n, where the i*'' element 
Pi denotes the transmit power of the i*'' transmitter The 
interference is treated as additive noise throughout this paper 
Ri denotes the maximum reliable rate of communication for 
the transmit-receive pair. Therefore, the achievable rate for 
the I*'' transmit-receive pair is written as: 



i?,(P) =l0g2 1 



(1) 



The next section finds the achievable rate region for such n 
transmit-receive pairs. 

III. Achievable Rate Region Frontiers for the 
Interference Channel 

First, the section considers the 2— user interference channel. 
The rate region problem is analyzed by formulating its under- 
lying nonconvex power control problem, for which we find a 
closed form analytical solution. The rate region is described by 
finding the maximum possible data rates achievable when each 
user is subject to a maximum transmit power constraint. This 
section then introduces the 3— user case to study the effect of 
adding a new dimension; and finally, by induction, the result 
is generalized for the ri— user case. 

A. 2-user Achievable Rate Region Frontiers 

In the case of the 2— user interference channel, Eq. ([T]l can be 
expressed as a function of Pi and P2 as P2), i = 1,2. 

For notational brevity, the channel gains are normalized by 
the noise variance, specifically: a = 51 i/cr^, b = 51,2/0"^, 
c = 92.2 /<^n, and d ~ 92.1 /(^n- Ri and R2 can therefore be 
written as: 

aPi 

(2) 

i?2(Pi,P2) = log2 (1 ' 



i?l(A,P2) =log2 1 



1 + dPi . 

For notational brevity, $(pi,p2) denotes a point in the rate 
region marked by having Pi = pi and P2 ~ P2- Effectively, 
the a;— coordinate of $(pi,p2) is Ri{pi,p2) = ri, and the 

coordinate of $(pi,p2) is -^2(^1,^2) ~ ^'2- The first 
objective is to find the achievable rate region frontiers of Eq. (|2]l 
through power control of Pi and P2, where each transmitter is 
subject to the maximum power constraint of Pmax- The frontier 
herein denotes the line (or generally, the (71 — 1) -dimensional 
surface for the 71— user channel) which traces the rate region 
via power control. 
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Fig. 2: 2— user power-control rate region 



B. Rate Region Frontiers Formulation 

The rate region frontier can be traced by setting Ri to a 
certain value ri, and then by sweeping ri over its full possible 
range from to i?i(-Pmax7 0) while finding the maximum R2 
value that can be achieved for each ri. From Eq. (|2]i, i?i is 
monotonically increasing in Pi and monotonically decreasing 
in P2, thus point i?i(Pmax, 0) corresponds to point $(P,„ax, 0) 
on the axis in Fig. |2] representing the maximum value Ri 
can attain. Similarly for the y— axis, the maximum value that 
i?2 can attain is i?2(0,Pmax), alternatively corresponding to 
point $(0, Pmax)- Those points represent the cases in which 
one of the users is silent, while the other is transmitting at full 
power Similarly, point $(Pmax, ^'max) has the coordinates of 

Rl{Pma^,Pma^) and R2{Pn 

Hence, for a constant rate i?i 



max; -fmax) 



ri, 



i?i(Fi,F2) =ri =log2 1 



aPi 



l + bP2 



(3) 



Therefore, the relation between Pi and P2 is obtained as 
follows: 



Pi = -(1 + &P2)(2'^^ -1). 



(4) 



From Eq. Q, for a constant i?i(Pi, P2) = n, i?2(Pi, P2) can 
now be written as a function of one parameter as i?2(P2), for 
Ri = ri, specifically: 



i?2(P2)=l0g2 1 



cP, 



l + -(l + feP2)(2''i -1) 



(5) 



It is important to analyze the behavior of i?2(P2) in terms of 
P2. This is presented in the following lemma: 

Lemma 1. Setting Ri at a constant rate, Ri[Pi, P2) = ri, 
P2(P2) is a monotonically increasing function in P2. 

Proof: The proof is provided in Appendix |A] ■ 



Using this lemma, the following corollary of uniqueness 
property is obtained: 

Corollary 1. For every rate tuple (?'*,?'2)> there is a unique 
power tuple ipi,P2)- 

Proof: i) a direct implication of monotonicity in Eq. (|5]l is 
that if R2 is equal to a constant at the rate of Pi = ri, then 
there is a unique P2 that achieves rj, ii) when P2 is determined, 
then Pi = pI is uniquely defined from Eq. (|4|, iii) from pi 
and P2, Ri is uniquely defined as Pi = ri = rl from Eq. (|3]l. 
Thus, pI and P2 uniquely define a point in the rate region with 
coordinates r* and rj. ■ 
In other words, any point in the rate region is achieved 
solely by a unique power tuple. This leads to what we denote 
by potential lines $ in the rate region, which are formed by 
holding one power dimension constant to a certain value and 
sweeping the other power dimension over its full range. In 
that regard, to describe a potential line marked by having Pi 
held at a constant power Pest, we use the following notation 
*^(-Pcst,:) to be equivalent to $(Pcst,P2) where Pi = Pest 
and < p2 < Pmax- Based on the uniqueness property 
just discussed, we have the following corollary (illustrated in 
FiglD: 

Corollary 2. Potential line^ along one power dimension do 
not intersect, i.e., $(:,p2) ond $(;,P2) do not intersect if P2 7^ 

P'2- 

The problem of finding the power control rate region fron- 
tiers then simplifies into finding the maximum value P2(P2) 
that can be achieved for any value of Pi(Pi,P2) = ri. 
Effectively, the formulation of the power control rate region 
is: 



arg max 

P2 



P2(P2) 



subject to Pi(Pi,P2) 
P < P 

^ 1 ^ -* max 



(6) 



1,2. 



ri is swept over the full range of Pi, i.e., < ri < 
Pi (Pmax, 0). The power control optimization problem in 
Eq. (|6]l is not straightforward as it is a nonconvex problem 1261 . 
However, by splitting the rate range of Pi into two intervals: 
Interval 7 for < ri < Pi (Pmax, Pmax) and Interval 2 for 
Pi (Pmax, Pmax) < < Pi (Pmax, 0), we are able to find a 
closed form analytical solution for each interval. The analysis 
of the optimization problem in Eq. (|6]l over these two intervals 
follows in the next two subsections. 



C. R2 Frontier for Interval 1: < ri < Pi (Pmax, Pmax) 

As Eq. (|3]l is monotonically increasing in Pi and monoton- 
ically decreasing in P2, ri can only exceed Pi (Pmax, Pmax) 
when P2 is less than Pmax- Thus, P2 = Pmax is attainable 
only when < ri < Pi (Pmax, Pmax), and P2 needs to 
be less than P,nax otherwise. From the proof provided for 
Lemma [T] where Eq. ^ is proved to be monotonically 
increasing in P2, and for the following Interval 1 range of ri. 

'The property in Corollary |2] is the reason for denoting these lines as 
potential lines, where the nomenclature is bonowed from electromagnetics 
based on a similar property for equipotential lines of an electric field 1251 . 
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Fig. 3: potential lines illustration in rate region and power 
region 



< ri < i?i(-Pmax, ^max), the solution to the optimization 
problem is: 



argmaxi?2(f2) = P^r 
Pi 



(7) 



Therefore, in this range of r\, using Eq. (01 and Eq. (|7]i, i?2 
is expressed by a function of r\ as follows; 



i?2(ri) =log2 1 



rP 



l + -(l+6P„,ax)(2'^^-l) 

a 



(8) 



Over this Interval 1 range of ri, the relation found in 
Eq. ([8]) describes the expression governing the potential line 
^(:,^inax), in which P2 is held at constant maximum power 
Pmax and Pi sweeps its full range from to Pmax- For brevity, 
the potential line $(:,Pmax) is denoted as $2, where the 
second power dimension, P2, is held at the maximum power. 
In this range of ri, $2 defined by Eq. ^ represents a power 
control frontier of the rate region as shown in Fig|2] 



p 

max ; ^ max 



) < ri < 



D. R2 Frontier for Interval 2: Ri(P„ 

Pi (Pmax, 0) 

Using symmetry of the previous result, for a constant rate 
P2 ~ ^2, there is a linear relation between Pi and P2. Thus, 
Pi(Pi,P2) can be written in function of one parameter Pi as 
follows: 



Pl(Pl)=l0g2 1 



aPi 



l + -(l + fiPi)(2'-^-l) 
c 



(9) 



By symmetry of the result in Lemma [T] Pi (Pi) is monoton- 
ically increasing in Pi. Thus, by symmetry, for the following 
range of r2: < r2 < P2(Pmax, Pmax), we have: 

argniaxPi(Pi) = Pmax- (10) 

Basically, the value found in Eq.(fTO]i describes the frontier for 
the following rate ranges of: < 7*2 < P2(^'max5 ^niax) and 
Pi (Pmax, Pmax) < 7^1 < ^i(J^max,0) - similar to the former 
subsection IIII-CI where P2 = Pmax described the frontier for 
the following rate ranges of: < ri < Pi (Pmax, Pmax) and 

P2(Pmax,Pmax) < ^2 < P2(Pmax,0). 

Therefore, the value of Pi at the frontier is: 

n P 



i?l(^max,^'2) =log2 1 



l + bP, 



(11) 



Hence, for this Interval 2 range of ri : Pi (Pmax, Pmax) < ^'i ^ 
Pi (Pmax, 0), the value of P2 that achieves the frontier follows 
as: 



P2 



aP„ 



2''i - 1 



1 



(12) 



So effectively, the value found in Eq. ( fT2b is the answer for the 
optimization problem in Eq. ^ for this range of R. Explicitly, 

for Pi (Pmax, Pmax) < ^-l < Pi (Pmax, 0), WC haVC: 



r> / r> \ ^ f ^-^max 

argniaxP2(P2) = - 1^1 



1 



1 



(13) 



Therefore, in this range of ri, using Eq. (fT2l i and Eq. (|5]l, P2 
is expressed in function of ri as follows: 

c 



P2(ri) =log2 1 



(aPm 



(2^^ - 1)) 



- 1) (1 + rfP„: 



(14) 



The relation found in Eq. (fT4l l describes the expression govern- 
ing the potential line $(Pmax, 0' where Pi is held at a constant 
maximum power Pmax and P2 sweeps its full range from to 
Pmax- Similarly for brevity, the potential line ^(PmaxiO is 
denoted as $1. In this range of ri, $1 as defined by Eq. (fl4l l 
represents a power control frontier of the rate region as shown 
in Fig|2] 

E. The 2— User Achievable Rate Region 

This subsection consolidates the two earlier results to fully 
describe the rate region frontiers. 

. For Interval 7 of ri: < ri < Pi (Pmax, Pmax) 



argmaxP2(P2) = Pa 
P2 
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Fig. 4; rate regions 7^l = 7^{$l} and 7^2 = 7^{<^>2} 



and the power control frontier, $2 
expressed as: 

/ 

i?2(ri) = log2 



$(:,-Pmax), is 



1 



rP 



\ l + -(l+6P„,ax)(2'- 

\ a 



(15) 



Let 7^.2 = Tl{^2} denotes the rate region outer-bounded 
by $2 as shown in Fig. |4] 

For Interval 2 of ri: i?i(Pmax, -Pmax) < < 

i?l(/'max,0) 



argmaxi?2(^2) = t 
P2 



1 / aP, 



2''i - 1 



and the power control frontier, $1 
expressed as: 

c 



R2in) = l0g2 1 



(2'-i -l)(l+dP„,ax) 




IS 



(16) 



Similarly, let TZi = 7^{$i} denotes the rate region outer- 
bounded by $1 as shown in Fig|4] 
In FiglH point A denotes point $(0, Pmax) (user 2 transmit- 
ting solely at full power), point B denotes point <l>(Pmax, Pmax) 
(both users are transmitting simultaneously at full power), and 
point C denotes point $(Pmax,0) (user 1 transmitting solely 
at full power). 

The rate region for a 2— user interference channel achieved 
through power control is obtained as: 



7^l u7^2. 



(17) 



Finally, the 2— user rate region, denoted as TZ, is found as the 
convex hull of the power control rate region. It is defined as: 



n = Convex Hull{7ei U 7e2}. 



(18) 



The treatment of the achievable rate region for the n— user 
interference channel follows next. It starts by considering a 
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Fig. 5: 3-user interference channel achievable rate region 

3— user interference channel to show the effect of adding a 
new dimension, and then generalizes the result for the n— user 
case. 

F. 3— User Example: Effect of Increasing P3 from to Pmax 
The rate region for the 3— user case is illustrated in 
FigS The following notation of $(^1,^2,^3) denotes a 
point in the rate region with coordinates of [i?i(Pi, P2, P3), 

P2(Pl,P2,P3), P3(Pl,P2,P3)]. Accordingly, $(:,Pmax,P3) 

describes a line characterized by sweeping the transmit power 
Pi of the first transmitter from to Pmax, with the second 
transmitter transmitting at Pmax and the third transmitter 
transmitting at a power value of P3. Similarly, Pi(:,Pmax, 
represents a surface in the rate region marked by sweeping the 
full range of Pi and P3, and holding P2 at P,„ax- 

When P3 = 0, the same setup and results that are described 
earlier in this section applies. Thus, for the rate range of < 
n < Ri{ 0) and < r2 < P2(0,Pmax,0) and 

P3 = 0, the frontier can be described as $(:,Pmax,0), which 
is the potential line from point A to point B in Fig. |5] As P3 
increases, the goal is to describe its effect and how it is traced 
in the rate region. 

Revisiting Eq. ([U, a fixed P3 has the effect of just an additive 
noise term in Pi(P) and P2(P)- Hence, all the previous 
results in section |lll] are applicable for any value of P3 in 
describing the frontier for Pi and P2, since the effect of 
P3 can be lumped in the noise term. Thus, for the range of 

< ri < Pi (Pmax, Pmax, Pa) and < r2 < P2(0,Pmax,P3), 

where P3 is constant, the frontier line on Pi and P2 is 
'!'(:, Pmax, P3), i-e., characterized by having P2 = Pmax- 
Consequently, the potential lines (or surfaces) concept in the 
3— user case carries through. 

Next, the frontier on P3 is described. For each value of 
P3, $(:, Pmax, P3) traces one of the highlighted curves in Fig. 
|5] For the collection of lines to form a frontier, we want to 
prove that at each increasing value of P3 these non-intersecting 
potential lines monotonically increase in the P3 dimension. 
This is evident from the relation between P3 and P3 in Eq. ([TJ. 
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The maximum value of R3 that can be achieved in this case 
is when P3 = P,nax, i.e., Rsi-, Pma^, Pmax)- Therefore, the 
highhghted frontier surface in Fig. |5] is the closed potential 
surface $(:,P,„ax,0- The boundary contours of this surface 
are the potential Unes: A^B,B^C,C^D, and D ^ A, 

defined as $(:,Pmax,0), $(Pmax, i^max, :), $(:, -Pmax, i^max), 

and $(0,Pmax, respectively. 

By symmetry of interchanging Pi, P2, and P3, the 3— user 
rate region is found via the convex hull onto the union of 
the regions bounded by these three surfaces; $(Pi„ax, 0' 
^{■, Pmax, 0' ™d $(:,:, Pmax)- The rate region 7?. is therefore 
expressed as: TZ = Convex HulljT^i U TZ2 U Ti-s}, where 
TZi ~ TZ{^i} is the region outer-bounded by the potential 
surface $i, where $i = $(..., P, = Pmax, • • ■) is the 
surface characterized by having Pmax in the i*'* power position. 
(Note that the intersection of potential surfaces is a potential 
line, as two of the dimensional inputs become equal, i.e., 

$(^'max,^'max,:) € 111 and $(Pmax,i^max,:) G 1^2-) 

G. n—User Generalization 

The case for n— user generalization is done by induction. 
For the n*'' added dimension to the existing 71 — 1 dimensions 
problem, the additional power effect of P„ can be lumped 
in the additive noise term of the existing expressions, and 
thus the results for Pi, ... , P„_i hold and carry through. The 
potential hyper-surfaces for fixed P„ are non-intersecting and 
monotonically increasing in P„, and thus the maximum outer 
limit is reached with P„ = Pmax for the appropriate range of 
i?i, . . . , Rn-i- Invoking symmetry, we can generalize over all 
the rates ranges, therefore arriving to the following theorem. 

Theorem 1. The achievable rate region of the n^-user inter- 
ference channel by treating the interference as noise is: 

n = Convex Hull{U^Li7?.,:}, (19) 

TZi = where $i is a hyper-surface frontier of n — \ 

dimensions, characterized by holding the i*'' transmitter at full 
power. 

Note that Theorem [T] also holds for different thermal noise 
levels or different maximum power levels. 

IV. Convexity Characteristics of the Power 
Control Frontiers For the 2-User Interference 
Channel 

This section focuses on the 2— user interference channel and 
studies the behavior of the power control frontiers, i.e., the 
potential lines $1 and $2, in terms of convexity and concavity 
in order to determine when the convex hull operation entails 
employing time-sharing. This happens whenever any of the 
potential lines, or segment thereof, is convex, which enables 
higher data rate to be achieved using time-sharing rather than 
using power control. Furthermore, specific results pertaining to 
the symmetric channel are presented at the end of this section. 

The power control frontiers equations for the 2— user case 
are: 



*i: R2{ri) = log2 1 



- (aPmax - (2^^ - 1)) 
(2'1 -!)(!+ dPmax) 



$2: R2{ri) log2 



1 



rP 



V 



1 



a 



bPrr 



0(2- 

It is not clear when $j is convex or concave, or whether 
it can exhibit a non-stationary inflection point. The non- 
stationary inflection point happens when the potential line has 
simultaneously a convex segment and a concave segment. The 
convexity behavior is thus treated next in more details. 

A. Convexity or Concavity of the Power Control Frontiers 

By using Eq. (|4]l when P2 = Pmax, the potential line $2 
depends on Pi through the following relation of ri and Pi: 

Pl = i(l+6Pmax)(2'^^-l). 

a 

Therefore, the second derivative of $2 with respect to ri leads 
to the following expression in function of Pi: 

= (a + adPiY - {a ~ a.){a — a + acPmax), (20) 



dri 

where a = d + dhP„ 



If the potential line is concave (i.e. ^ < 0) then the 
enclosed region Tl{^2} is convex. The rate region is defined 
to be convex when a straight line connecting any two points 
inside the rate region is entirely enclosed in the rate region. In 
contrast, if the potential line is not concave, i.e., if it is convex 
or exhibits a non-stationary inflection point, then we describe 
its enclosed region as being concave; as in this situation, the 
aforementioned definition of a convex region does not hold. In 
summary, if is concave, then 7^{$i} is convex, and 7?,{$i} 
is concave otherwise. 

Let 5R(-) be the real operation, the inflection threshold Qi is 
defined as: 



(Q ~ - ce + QcPmax)) - a 
ad ' 



(21) 



where Qi was derived such that sign ^^^^ = 
sign(Pi — Qi). Therefore, it suffices to study convexity 
or concavity of potential line $2 by examining the sign 
of (Pi — Qi). Note that the inflection threshold Qi only 
depends on system parameters a, b, c, d, and Pmax- The 
relation is nonlinear By plugging in the respective values, the 
convexity behavior is assessed, and it can be decided whether 
time-sharing is needed. 

Thus, $2 can exhibit the following convexity behaviors: 

• Qi > Pmax- then Pi — Qi < for all the range of Pi, 
thus $2 is concave. Operating via power control is optimal 
in leading the highest achievable data rate, and no time- 
sharing is needed. See $2 in case (i) in Fig|6l 

• < Qi < Pmax: $2 exhibits a non-stationary inflection 
point when Pi = Qi (see point D in Fig|6]and Fig|7]i. In 
this case: 

- for < Pi < Qi: line $(0 ; QijPmax) is concave, 
i.e., the potential line segment ^ad is concave as in 
Fig|6] case (ii) and Fig|2] (a). Operating via power 
control to trace the segment $(0 : Qi,Pmax) is 
optimal. 
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Fig. 6: Pmax = 1; (i): concave $2 and $1 {[a,b;c,d] ~ 
[10, 1; 1, 10]), (ii): $2 with inflection point D, concave $1 
([10, 1; 4, 10]), (iii): convex $2, concave $1 ([10, 1; 6, 10]), (iv): 
convex $2 and $1 ([10, 15; 4, 10]) 



- for Ql < Pi < Pmax: line $(Qi : Pmax,Pmax) 

is convex, i.e., the potential line segment is 
convex as in Fig|6]case (ii) and Fig|7](a). Therefore, 
operating via the time-sharing segment between the 
inflection point D and point B is optimal. 

. Qi < 0: then Pi - Qi > for all the range of Pi, 
thus $2 is convex. See <I>2 in cases (iii) and (iv) in Fig|6] 
Depending on the a, b, c, d, and Pmax parameters, it is 
optimal to apply time-sharing with the following options: 

- between point A and point B, see Fig]?] (c). 

- between point A and a point on the concave segment 
of $1, see Fig|2](b). 

- between point A and point C, see Fig|7] (d). This 
is effectively Time Division Multiplexing (TDM), 
where each user transmits solely at any point of time. 
This form of dimension-orthogonality occurs when 
the interference is very strong rendering the cost 
too high for having simultaneous transmission. The 
upcoming subsection II V-BI explores the optimality of 
operating via TDM, or equivalently, it explores when 
this cost is deemed too high. 

By symmetry, the potential line $1 exhibits similar convexity 
behavior: a) it is convex when Q2 < 0, b) it is concave when 
Q2 > Pmax, and c) it exhibits a non-stationary inflection point 
when P2 = Q2- Hereby, with /? ~ {b + bdPynax), the inflection 
threshold Q2 is defined as: 



^(Vic - P){C - /? + QcPnax)") - 

cb 



(22) 



Note that by virtue of how the inflection threshold is situated 
with respect to power P;, whenever any of the potential lines 
exhibit an inflection point, the order of which segment is 
convex or concave is not arbitrary. It always starts concave 
in the segment closer to the coordinate axis. Explicitly, if $2 
and $1 exhibit inflection points, then tracing the rate region 




(b) 



A 


(d) 











Fig. 7: types of time-sharing. 



from left to right: $2 is concave then transitions to convexity, 
and $1 is convex then transitions to concavity. 

B. When is TDM Optimal? 

Discounting the case when $1 or (E>2 exhibit non-stationary 
inflection point for simplicity, and focusing on the case when 
both potential lines are convex (i.e., Qi < and Q2 < ), it 
is important to know when TDM is optimal. Under the afore- 
mentioned assumption, this translates to determine when time- 
sharing between point A and point C is better than time-sharing 
through intermediate point B. This is done by comparing the y- 
axis ordinate of point B, P2(Pmax, Pmax), relative to the y-axis 
ordinate from the straight line connecting points A and C at 

ri = Pi (Pmax, Pmax), denoted as PP(ri)|ri=_Ri(P,„ax,Pmax)' 

as shown in Fig]?] (d). Namely, TDM is optimal when 

^P(^l)lri=Pi(P,„,„,P,„,x) > P2( Pmax, Pmax); that is: 

-log2(l + cP„,^) 



log2(l+aP„ 



■log2(l 



aP„ 
1+bPm 



+ l0g2(l +cP,nax) 



This leads to the following Lemma |2] 

Lemma 2. Operating via TDM (i.e. one transmitter solely 
transmitting at a certain time) is optimal in achieving the rate 
region when 



(l+cPmax)(l+rfP„ 



> 



bP 

max 



(23) 



with 7 = log2(l + cPmax)/ log2(l + aP,„ax)- 

Note that the condition found in Eq. i23[ is a nonlinear 
relation between the interference channel variables of a, b, c, d, 
and Pmax- This motivates the following subsection II V-C I to treat 
the 2— user symmetrical channel. 

C. Symmetric 2— User Interference Channel 

This subsection treats the symmetric 2— user channel, mainly 
analyzing the expression in Eq. ( |23] | in order to derive clear 
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Fig. 8: symmetric 2— user rate region: Pmax = 1, a = 1, b^'^ = 
0.2, = 0.4, = 1, = V^, and fot") = 3. The 

threshold b* to switch to time-sharing is equal to \/2 from 
Eq. (El. 



insights. For the symmetric channel, a = c and b = d, the 
expression in Eq. ( |23] | simplifies, and leads to the following 
corollary on the TDM optimality condition: 

Corollary 3. For the symmetric 2— t/ier interference channel, 
operating via TDM is optimal in achieving the rate region when 



b> 



p 

^ max 



(24) 



In other words, when the interference is weak (i.e., b is 
below the threshold in Eq. (l24l i). then it is best for both 
transmitters to transmit at full power When the interference 
increases and exceeds the threshold in Eq. ( |24| ). then the TDM 
scheme becomes optimal. In this scenario, the users can no 
longer share the same resource, and thus they have to use it in 
an orthogonal fashion. An example is illustrated in Fig|8] The 
application of the usage of the different types of time-sharing 
and the advantage over power control has been discussed in 
ll27l . when it was applied to a 2— sector interference channel 
in a cellular setting. 

Remark: For high SNR (i.e. aPmax S> 1), Eq. (l24l) reduces 

> 



to &Pmax > -y/aPmax, which coincidcs with the results in 
where this condition marks the interference power threshold 
above which treating the interference as noise is no longer 
optimal in the Degrees of Freedom sense and no longer within 
a gap of 1-bit from achieving capacity. 

In addition. Appendix |B] proves that the expression in 
Eq. (l24l) is a sufficient condition for both frontiers $i and 
4>2 to be convex, i.e., Qi and Q2 are always < 0; which is the 
starting necessary condition of subsection IIV-BI when treating 
the general asymmetric channel. Therefore, when b satisfies the 
expression in Eq. (l24l i. the frontiers potential lines are always 
convex and TDM is optimal. 

The extension of Corollary [5] to the 71— user symmetric 
interference channel is provided in Appendix |Cl leading to 




Fig. 9: 2— user crystallized rate region: (a) rate region through 
power control, (b) crystallized hull overlaid on top of the rate 
region 



the following TDM optimality condition when: 



b > 



aP„ 



(1 + aPmax) 



1/r. 



- 1 



1 



(n - 1)P„ 



(25) 



V. Crystallized Rate Region 

Based on the discussion in the previous two sections, this 
section introduces a novel approach into simplifying the rate 
region in the n— dimensional space by having only an On/Off 
binary power control. This consequently leads to 2" — 1 corner 
points within the rate region. By forming a convex hull through 
time-sharing between those corner points, it thereby leads to 
what we denote a crystallized rate region ||2]. The concept 
of the crystallized rate region has since been extended to 
the MIMO and the MIMO-OFDM channels El, ||29l. This 
section focuses on the crystalUzed rate region formulation and 
its evaluation. 

As illustrated in Fig|9]for the 2— user case, the crystallized 
rate region is an approximation of the achievable rate region 
formed by a convex hull of straight lines connecting points 
A, B, and C. Those corner points are achieved through binary 
power control with the transmitters employing either zero or 
full power. For the 2— user interference channel, there exist 3 
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g=[3 1 0.5; 1 2 0.5; 0.5 1 2], P =1 

^ ' ' ■" max 




g=[3 1 0.5; 1 2 0.5; 0.5 1 2], P =1 

^ ^ max 




Fig. 10: 3-user crystallized rate region: (a) time-sharing crys- 
tallized hull, (b) crystallized hull overlaid on top of the power 
control rate region 

corner points; similarly for the n— user case, there exist 2" — 
1 corner points. Note that in Fig|9] the time-sharing convex 
hull connects point A and point C though point B; whereas if 
the interference is strong, the convex hull is formed by time- 
sharing A and C only. Fig[TO] shows an example of the 3— user 
case where Fig{TO](a) is the crystallized hull, and Fig{TO](b) is 
the crystallized hull overlaid on top of the power control rate 
region. In the 2— user dimension, the time-sharing convex hull 
is a set of straight lines connecting two points. In the 3— user 
dimension, it is a set of polygon surfaces connecting three 
points, see Fig. [TT] In the n— user dimension, each polygon is 
formed by connecting n points, and hence the number of such 
polygons is ^^). 

A. System Time-sharing Coefficients and New Rates Equations 

Recall that with power control, the achievable rates for the 
2— user interference channel are given in Eq. (|2]i. The paradigm 
with the crystallized rate region approach is different. Instead 
of formulating the problem as a power control problem for 
finding Pi, the crystallized rate region formulation focuses 
on finding the appropriate time-sharing coefficients of the 



g=[3 1 0.5; 1 2 0.5; 0.5 1 2], P =1 




Fig. 11: a desired operating point in a time-sharing polygon 

2" — 1 corner points. For the 2— user case, let 9~ [9i, 02, ^'a]"^, 
J2i Si — 1, denote the system time-sharing coefficients vector 
of respective corner points $(Pmax,0) (user 1 transmitting 
solely with a time-sharing coefficient 6i), $(0, Pmax) (user 
2 transmitting solely with a time-sharing coefficient 62), and 
*I'(-fmax, -Pmax) (both users transmitting with a time-sharing 
coefficient 6*3). The reason for labeling a system time-sharing 
coefficients vector is to emphasize the combinatorial element 
in constructing the corner points, where the cardinality of 
|6»| = 2" - 1. 

For 2— user case, in contrast with Eq. the rates Ri and 
i?2 for the crystallized region are: 

Ri{e) = 9i log2(l + aP„,ax) + O3 log2 (1 + I'fbpL. ) ' 

i?2(0) = 92 l0g2(l + cP„,ax) + ^3 loga (l + T^^T^) ' 

""^26) 

Any solution point on the crystallized hull lies somewhere on 
a time-sharing line connecting two points for the 2— user case; 
and similarly in the 3— user case, the solution point lies some- 
where on a time-sharing plane connecting three points (see 
FigfTTTl. then by induction we obtain the following corollary: 

Corollary 4. For any solution point on the n~user crystallized 
rate region, the system time-sharing vector 6 has at maximum 
n nonzero coefficients out of its 2" — 1 elements. 

Corollary |4] could also have been reached invoking the 
Caratheodory theorem 1301 . where, in the paper's context, any 
point lying on the hyper-surface of dimension n — 1 is enclosed 
in a convex set of n or fewer points inside the hyper-surface. 
Therefore, if those points are the corner points, any solution 
point on the n — 1-dimensional hyper-surface will be the result 
of time-sharing at maximum n corner points. 

Let a'*^' denote the transmitters power mask that character- 
izes the k*^ corner point, specifically: 

j = 1 , • • • , n, where i is the transmitter index, 
A: = 1, • • • , 2" — 1, where k is the corner point index. 

(27) 
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af^ is the action taken by the i*'^ user in characterizing the k*'^ 
corner point; the binary transmit action is either being silent 
or transmitting at full power. Let «_/ indicate the interferers 
transmit power mask, which is derived from vector a^^'> by 
excluding the i*^ transmitter's action. The generalization for 
the n— user crystalUzed rates equations now follows as; 



(fe) D 



k=l 



(28) 



where gu is the desired channel gain from the i transmitter to 
the i*'' receiver normalized by the noise power; gi is the noise- 
normalized channel gains vector received at the i*'' receiver 
from all the transmitters excluding the i*'^ transmitter The 
length of each vector a^^' and gi is n — 1. 



B. Evaluation of Crystallization 

In this section, we compare the crystallized rate region 
and the rate region bounded by the power control potential 
lines. In effect, we are evaluating how much gain or loss 
results from completely replacing the traditional power control 
scheme (see Eq. (|2]i) with the time-sharing scheme between 
the corner points (see Eq. (|26]|). For this purpose, we consider 
the symmetric channel, where a = l,Pmax = 1, and we 
increase the interference b to vary the signal to interference 
ratio SIR = a/b. Two metrics are used as a measure: (a) 
the area of the rate region, and (b) the maximum bit rate gap 
between the traditional power control and the crystallized rate 
region scheme. 

The values of the area bounded by the power control poten- 
tial lines and the area bounded by the crystallized rate region 
are plotted in Fig fT2l a). For weak interference, or equivalently 
noise-limited regime, point B is used in constructing the 
crystallized region. As the interference increases beyond a 
certain threshold level, time-sharing through point B becomes 
suboptimal, and time-sharing A-C becomes optimal. The exact 
switching point between power control to time-sharing is given 
in Eq.(l24li. In Fig fTSl a). this happens at the intersection of the 
time-sharing line through point B and the A-C time-sharing 
line. As indicated in Fig fT2l a). there is no significant loss in 
the rate region area if time-sharing is used universally instead 
of traditional power control; in fact in some cases time-sharing 
offers considerable gain. Specifically, whenever the potential 
lines exhibit concavity, time-sharing loses to power control; 
whenever the potential lines exhibit convexity, time-sharing 
wins over power control. Different values of a also lead to 
the same conclusion. 

The second measure, the maximum rate gap percentage 
between the rate achieved by traditional power control and the 
rate achieved by using the crystallized rate region is plotted in 
Fig fT2l b). The maximum bit rate loss of using the crystallized 
hull compared to power control does not exceed 1%, and the 
crystallized strategy is therefore quite attractive. It is arguable 
that as the interference becomes the network bottleneck there 
is little reason to implement non-binary power control, and that 
dimension orthogonalization becomes the primary objective. 
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Fig. 12: (a) comparison of the areas bounded by crystallized 
hull and power control, (b) maximum rate gap percentage of the 
difference between crystallized hull and power control rates, 
with "b" from -20dB to OdB 

VI. Conclusion 

The achievable rate region for the n— user interference 
channel was presented when the interference is treated as 
noise. The results were first found for the 2— user interference 
channel, then they were extended to the 3— user case to show 
the effect of adding an additional dimension. Subsequently, 
they were generalized for the n— user case. The n— user rate 
region was found to be the convex hull of the union of the 
n rate regions, where each rate region is upperbounded by a 
(n— 1) -dimensional hyper-surface characterized by having one 
of the transmitters transmitting at full power For the 2— user 
interference channel, those hyper-surfaces become what the 
paper refers to as power control potential lines. Hence, to 
determine when time-sharing should be used in forming the 
achievable 2— user rate region, the paper studied the potential 
lines convexity behavior. Finally, the novel concept of the 
crystallized rate region was introduced and evaluated. The 
crystallized rate region is described by composing a time- 
sharing convex hull onto the 2" — 1 corner points obtained 
from On/Off binary power control. The evaluation of the 
crystallized rate region shows little value in the implementation 
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of non-binary power control compared to the use of binary 
power control in conjunction with time-sharing. In effect, 
the crystallized rate region approach offers a new perspective 
of looking at the achievable rate region of the interference 
channel. 

Appendix A 

Proof that R2{P2) is Monotonically Increasing in 

P2 

Proof: Effectively, Eq. (|5]) is in the form of /(I + g{x)). 
As /(•) is monotonically increasing in its argument, it suffices 
to prove that g{x) is monotonically increasing in x. Therefore, 
define .9(^2) as: 

acP2 



9{P2) = 



dg{P2) 
dP2 



a + d{l + bP2){T^ - 1)' 
ac{a + d(2'~i - 1)) 



(29) 



(30) 



^2 (a + d(l + 5P2)(2'-i - 

The numerator in Eq. ( l30l l is nonzero if a 7^ and c 7^ 
(a = or c = are the trivial cases where the rate region is 
either a line or the point zero). As ri > 0, then {2^^ — 1) > 
0. Thus dg{P2)/dP2 is always > for non-trivial cases of 
a and c. Thus, g{P2) is monotonically increasing in P2, and 
equivalently i?2(^2) is monotonically increasing in P2- ■ 

Appendix B 

Proof that Eq. (1241) is a Sufficient Condition for 
Both $1 and $2 to be Convex 

Proof: For the symmetric case, Qi = Q2 = Qsym can be 
written as 



„ _ SR(V(a-^)(a-g + a^P.uax))-g 

^ sym — , 1 

ab 

where 9 = b + bi^ P^ax- Qsym can also be written in this form: 



Q 



sym 



ab ' 



where Ti = a — 9 = a — b — b^ Pj^ax, and T2 = a—O + a^P^ax- 
From the expression in Eq. (|24] |. a can be alternatively upper- 
bounded as a < {b^Pmax — 1/^'max)- Therefore, Ti is upper- 
bounded as: 



Tl < -1/Pn- 



(31) 



From Eq. ( [3T1 ), Ti is always negative. T2 however can be 
positive or negative, evaluated as follows: 

• 72 > 0: then 3?(VTi72) = 0, and as 9 is always positive, 

then Qsym < 0- 
« r2 < 0: 5R(\/7iT^) is > 0. In this case, the numerator of 
Qsym can be written as: 

num 



(Qsym) = Vi^-aW-a-a^Paax) - 9. 

Given the fact that (6* — a — a?Praax) < (9 — a), then 
num{Q sym) can be upper-bounded as: 



num{Qsym) < -9<-a<0. 
Hence, the frontiers potential lines $i and $2 are convex. 



Appendix C 

Extension of Corollary [3]to the n-usER case 

Focusing on the (_Ri,i?2) 2D rate region of an n— user 
symmetrical rate region, if Eq. (l24l i is satisfied. Appendix |B] 
proved that it is a sufficient condition to make sure that $1 
and $2 are convex and that they do not have an inflection 
point. Let 63 denote the threshold value in Eq. (|24] |. Hence, if 
the interference channel gain b > then $1 and $2 in the 
(i?i,i?2) 2D rate region are convex. As more users transmit, 
they will cause additional interference power on i?i and i?2. 
Taking the 3— user case as an example, if the 3'"'' user transmits, 
then i?i(Pniax,-P2,-P3) = log2(l + T+febft)- Projecting 
the 3D rate region into the (i?i,i?2) 2D plane will lead to 

Rl{Pn.ax,P2) = log2(l + ^^), whcrC b' = b+bPz/P2, 

which is greater than b. Thus, if the (i?i,i?2) 2D rate region 
have a convex $1 and $2, then > ^2' and projecting higher 
71— dimensional rate region frontiers into the R2) 2D plane 
will always result in power control frontiers, denoted as 
and that are always convex - due to the fact that they 
can be obtained using an interference power gain b' where 
b'>b> b*. 

The goal of this appendix section is to find the n— user 
symmetrical interference channel gain threshold, denoted as 
6* , such that when 6 > 6* , TDM is optimal in leading higher 
achievable rate region. Thus, as a starting assumption (which 
is proved to hold later on), assume 6* > 62,V7i > 2; which 
therefore leads to the property that all the ?i— dimensional 
hypersurfaces are convex - as the projections on all the pair- 
wise 2D planes result in convex power control frontiers. 

With a TDM solution, users transmissions are orthogonal in 
time, and therefore, whenever user i transmits all other users 
are silent. The maximum TDM rate achievable for user i is 
when Pi = Pmax, and in the context of an n— user symmetrical 
channel, it is equal to Rj = R^^^^ = log2(l + aP,nax), and 
Rj = 0,Vj 7^ i. Therefore, whenever user i transmits under 
TDM, its maximum achievable rate in the n— dimensional rate 
region is a point on the z*'* axis with a value equals to i?^^^^; 
let ^f^^^ denotes such a point. 

Let "H be a hyperplane formed by connecting via time- 



sharing the $ 



TDM 



1, • • • ,n points. Define the origin 



point O as the point with coordinates Rf — 0,Vi. Define 
point B on the power control achievable rate region when 
all users transmits at the same time; thus point B coordinates 

areiif ^1 ■ 



log2 



, Vi. Define point B' 



\^ ^ 1 + {n - l)bP„ 
on the time-sharing hyperplane H when all users transmits in 
TDM for an equal amount of time. Let OB denotes a vector 
from point O to point B, and similarly, OB' denotes a vector 
from point O to point B'. Let || • || denotes the vector length; 
for instance, for OB, ||OB|| is the distance from point O to 
point B. Therefore, a time-sharing TDM solution is optimal 
when 1 1 OB' 1 1 > ||OB||; which has the interpretation of the 
hyperplane H leading higher achievable rate region than the 
power control hyper-surfaces frontiers. 

Based on the aforementioned coordinates of point B, 
||OB|| = x^log2(l + ,+^:'^T)bP^J . For IIOB'II, it is flie 
shortest distance from point O to the hyperplane H; which is 
equal to the absolute value of the dot product of the unit normal 



12 



vector of % and the vector OK, where point K is a point on 
%. The unit normal vector of % is equal to -^[1, • ■ ■ , 1]"^- 
By choosing point K to be equal to the point on the 

i?i axis, ||0$™*^|| = [i?^^^^ 0, • • • , 0]^, where R^^^' has 
been defined earlier to be equal to log2(l + aPniax)- Therefore, 

IIOB'II = ^l0g2(l + aP,„ax). 

Hence, a time-sharing TDM solution is optimal when 



^ log2(l + aP„, 



> 1 



aP„ 



1 + (n - l)bP^ 



(32) 



leads to Eq. (EU. The threshold 6; 



- 1 j („_]^|p — is monotonically increasing 
which starts at Eq. (l24l i for n — 2 and 



Expanding Eq. 

/ aP,„ax 

V(l+a-Pmax)l''"-l 

for a positive n, 
flattens out asymptotically for large n. Using the first order 
approximation of (1 + aP^ax)^/" to 1 + ^ ln(l + aPmax), 
asymptotically converge to: 

lim &* a/ ln(l + aPmax), 

independent of the number of users n. 
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